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We study the transverse and longitudinal momentum spectra of fermions produced in a strong, 
time- dependent non-Abelian SU(2) field. Different time-dependent field strengths are introduced. 
The momentum spectra are calculated for the produced fermion pairs in a kinetic model. The 
obtained spectra are similar to the Abelian case, and they display exponential or polynomial be- 
haviour at high pt, depending on the given time dependence. We investigated different color initial 
conditions and discuss the recognized scaling properties for both Abelian and SU(2) cases. 
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I. INTRODUCTION 

During last years large amount of data on particle 
spectra have been collected in relativistic heavy ion col- 
lisions at the Super-Proton Synchrotron (SPS, CERN) 
and at the Relativistic Heavy Ion Collider (RHIC, BNL) 
at cm. energy of = 10 - 200 GeV [i, 2,] in a 
wide transverse momentum range, < pt ^ 20 GeV. 
Since the microscopic mechanisms of hadron production 
in hadron-hadron and heavy ion collisions are not fully 
understood, thus it is very important to improve our the- 
oretical understanding on this field. The forthcoming 
hadron and heavy ion experiments at the Large Hadron 
Colhdcr (LHC, CERN) at = 5500 GeV will increase 
the transverse momentum window to < pt < 50 GeV. 
Thus LHC experiments will become a decisive test be- 
tween different perturbative and non-pcrturbative mod- 
els of hadron formation, especially in the high-py region. 

Theoretical descriptions of particle production in high 
energy pp collisions are based on the introduction of chro- 
moelectric flux tube ('string') models, where these tubes 
are connecting the quark and diquark constituents of the 
colliding protons [1, 0, IB Hi- String picture is a good 
example of how to convert the kinetic energy of a col- 
lision into field energy. New hadrons will be produced 
via quark-antiquark and diquark-antidiquark pair pro- 
duction from the field energy, namely from the unstable 
fiux tubes. These models can describe experimental data 
very successfully at small pT^ especially at pr < 2 — 3 
GeV. At higher pr one can apply perturbative QCD- 
based models 0, B Q , which can provide the necessary 
precision to analyse nuclear effects in the nuclear colli- 
sions. 

However, at RHIC and LHC energies the string den- 
sity is expected to be so large that a strong collective 
gluon field will be formed in the whole available trans- 
verse volume. Furthermore, the gluon number will be 



so high that a classical gluon field as the expectation 
value of the quantum field can be considered and in- 
vestigated in the reaction volume. The properties of 
such non-Abelian classical fields and details of gluon 
production were studied very intensively during the last 
years, especially asymptotic solutions (for summaries see 
Refs. lOj llj)- Fermion production was calculated re- 
cently Lattice calculations were performed to de- 
scribe strong classical fields under finite space-time con- 
ditions in the very early stage of heavy ion collision, es- 
pecially immediately after the overlap [3, [3, [3] ■ New 
methods have been developed to investigate the influence 
of inhomogeneity on particle production [TtI . [TsI . [19| . 

Fermion pair production together with boson pair pro- 
duction were investigated by different kinetic models of 
pa.rticle production from strong Abelian [20l. [2ll . [2^ . [23l . 
lii, HE Hi, [13 and non-Abelian [H, [HI, infields. These 
calculations concentrated mostly on the bulk properties 
of the gluon and quark matter, the time evolution of the 
system, the time dependence of energy and particle num- 
ber densities, and the appearance of fast thermalization. 

Our main interest is the transverse momentum distri- 
bution of produced fermions and bosons, however the lon- 
gitudinal momentum can become equally interesting con- 
sidering wide rapidity ranges. In our previous paper (see 
Ref. ,33i] ) we investigated the Abelian case, namely parti- 
cle pair-production in a strong external electric field. We 
have focused on the numerical solution of a kinetic model 
and discussed the influence of the fleld strength apply- 
ing different time dependences. We have demonstrated 
a scaling behaviour in time and transverse momenta. 



namely t ■ Eq^^ and kx/E^^'^ . The kinetic equation and 
the numerical calculation yielded a fermion dominance in 
the mid-rapidity region. In case of realistic Bjorken type 
time evolution our numerical result on fermion spectra 
has overlapped the boson spectra obtained in 1-1-2 dimen- 
sional lattice calculations both in magnitude and shape. 
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In this paper we solve the kinetic model in the presence 
of an SU(2) non-Abehan color field. We focus on the de- 
termination of the appropriate kinetic equation system 
and its numerical solution in case of different initial con- 
ditions. Section 2 summarizes the kinetic equation for 
color Wigner function. In Section 3 we introduce phe- 
nomenologically the fermion distribution function. In 
Section 4 the kinetic equation and its simplifications are 
displayed in detail, especially with zero fermion mass. In 
Section 5 we present kinetic equation for the case, when 
external non-Abelian field have fixed direction in color 
space. The appropriate form of the kinetic equation is 
determined, which was solved numerically at different 
time-dependent strong fields. Results on field compo- 
nents and fermion distribution functions are displayed 
in Section 6. We display recognized scaling properties 
of the solutions in Section 7. In Appendix A we show 
a special non-Abelian setup, which leads to the similar 
results obtained in the Abelian case. 



II. THE KINETIC EQUATION FOR THE 
WIGNER FUNCTION 

The fermion production in a strong external field can 
be characterized by a space homogeneous Wigner func- 
tion, M^(k;t). The evolution of this Wigner function is 
investigated by the kinetic equation in the frame of the 
covariant single-time formalism, where a time-dependent 
Abehan (A^) [34, 35, 36, 37] or non-Abehan (A^) 28, lH 
external field is included. Here we choose a longitudinally 
dominant color vector field in Hamilton gauge described 
by the 4-potential 



Al = iO,~A-) = {0,0,0, A-^) 



(1) 



The Wigner function depends on the 3-momenta k = 
{ki, k2, k^). In the "instant" frame of reference ^2^ we 
obtain the kinetic equation for W^(k; t) as 

9tW+^—{A{W,Fo,}+ 
8 oki 

+2{F,„[M/,7V]}- [i^..,{W^,7V}]) = 
= ih{-f°-f\ W} - im[-f°, W] + ig [A, , [7%', W]] . (2) 

Here m denotes the current mass of the fermion produced 
in the strong field and g is the coupling constant. 

The color decomposition with SU{Nc) generators in 
fundamental representation [t"-) is 



W^W + We, a = 1,2,. 



1 



(3) 



where is the color singlet and is the color multi- 
plet component (triplet in SU{2) with Nc = 2). 
The spinor decomposition is the following: 

The time evolution of the Wigner function is described 
by eq. Using the color and spinor decomposition 



terms in eqs.Q-Q, we will determine the time evolution 
numerically. However, our main goal is to obtain the 
transverse momentum spectra of the produced fermion 
pairs, thus we need to introduce the fermion distribution 
function. There is no straightforward way to define this 
distribution function, but we can find a phenomenological 
way, summarized in the next chapter. 



HI. THE FERMION DISTRIBUTION 
FUNCTION FROM THE ENERGY DENSITY 

The energy density carried by the produced fermions 
is defined throughout the Wigner function [2^, [35[ : 



£f = Ti{{m-Yh)W + uj{k)) 



(5) 



where the second term fixes the zero energy level to the 
physical value (note, that expression ([5]) is free from di- 
vergent vacuum contribution). Here the one particle en- 
ergy is denoted by CLi(k) = \/l?~fn?. The trace is over 
color and spinor and the averaging is taken over mo- 
menta: 



(pk 



-A'(k;0. 



(6) 



After the color and spinor decomposition one obtains 
Ef = 2Nc{2ma' + 2\^h'' +UJ{\^)). 



(7) 



In parallel, focusing on one-particle distribution function 
the usual energy density formula has the following form: 



e; = 4iV,(^(k)/;(k) 



(Pk 

(2^ 



;(k)//(k). (8) 



Note, that both averages ([5]) and ([8]) may contain ultra- 
violet divergence, that can be properly regularized (for 
Abehan fields see e.g. [2^, HO]). However, since we are 
not going to calculate bulk properties of produced pairs, 
the issue of regularization is not important in current 
consideration and will be considered elsewhere. 

Now, combining eq.([7]) and eq.® one can derive phe- 
nomenologically a distribution function, namely 



//(k,t) 



ma''(k,t) -|-kb''(k,t) 1 
J(k) ^ 2 



(9) 



The time dependence in ff (k, t) is connected to the time- 
evolution of a^(k, t) and b''(k, t), which is followed by 
solving the decomposed kinetic equation of eq. ([2]). 

The fermion distribution function should be zero in 
vacuum. The expression in eq. ([9]) satisfies this request. 
Indeed the vacuum solution (A^ = 0) for the singlet 
Wigner function has the following form (see Ref. j28| ) 



1 771 

2 



- k7 



c^(k) 



(10) 



3 



The vacuum solution for the multiplet Wigner function is 
zero, = 0. The substitution of this Wigner function 
into the definition of the distribution function in eq.® 
leads to // = in vacuum, which is physically correct. 
Furthermore, one can see that //(k, t) is positive definite 
in presence of a non-zero field (A^ ^ 0). 

Although our fermion distribution function //(k, t) has 
been postulated phenomenologically and possibly it is 
not a unique solution, but it seems to be the right tool 
to investigate the time evolution of the energy distribu- 
tion at the microscopical level in a space homogeneous 
environment. 



IV. KINETIC EQUATION FOR THE WIGNER 
FUNCTION IN SU(2) 

Our basic aim is to determine the fermionic distribu- 
tion function and its time evolution in the presence of a 
strong, time-dependent non-Abelian SU(2) field. For this 
task we need to substitute the color and spinor decom- 
posed Wigner function of eq. ^ into the kinetic equation 
in eq. ([2]). After evaluating the color and spinor inde- 
ces we obtain a system of coupled differential equations, 
which consists of 32 components in SU(2): 



4 dk'i 

dta'^ + gE'^^a' = -4kc?, 
oh 



d 



oh 

dth" + -E^'—h" = 2k X -I- 4mcf , 
4 Oh 

dth" + gE^—h" = 2k x d° + 4mc? 
oh 

+2gr^''A%, 

dtcl + ^E-—c1 = ka^-mb^ 
4 dh 

dtcl+gE'^^cl - ka''-mb'' + 
oh 



(11) 
(12) 

(13) 
(14) 
(15) 

(16) 
(17) 



2gfabc^A>' X c^), (18) 



9 j^a 9 „ 



^ + 4 dh^ ~ ' 



(19) 



5tc^ + gE-—c^, = ke'^ - 2gr-{K'' x c^), (20) 
oh 



(21) 



4 Oh 

dtd^o+gE'^-^^d'o = 2me' + 2gr'''A''d', (22) 



dtd' + ^E''—d'' = 2kxb^ 
4 oh 



(23) 



dtd- 



d 

oh 



2k X b'^ + 2gt^''A°d\ 



9^a 



d 
dh 



-4kc? - 2mdn 



4 dh 



-4kc^ 



2mdg. 



0, (24) 
(25) 
(26) 



Here we use the notation 6^ — (6q , — b*^'"), and a simi- 
lar one for d^^"" . The antisymmetric tensor com pon ent of 

1. The 



_ Aa 



the Wigner function is defined as c^" = (ci, C2) 
strength of the color electric field is E'^ = E^ - 

Since the mass of light quarks is small [ill, we neglect 
the mass term in our SU(2) calculation. In this case 
the distribution function for massless fermions is solely 
defined by b'*: 

. ^ kb* 1 , , 

//(k.^) = ^+o' (27) 



The zero fermion mass leads to the simplification of 
the above kinetic equation, which finally is splitted into 
two independent parts: one for a*'°, c^'°, e**'" and 
another one for 6g °, b**'", dp °, d**'". One can recognize 
that the second set completely defines the evolution of 
the fermion distribution function. This way the chiral 
symmetry of the massless case is manifested itself. Thus 
we focus on these equations, only: 



4 Oh 



dtd^^ 
dtb^o 
dtd^o 



d 

-gE^^h^ 
oh 

lE'^^d" 



4 
gE'^ 



dh 
d 

dh' 



2k X d", 
2k X d° 
2k X b", 
2k X b" 



9 pa d 

4^ dh^" 

9 pa d 

4^ W/' 



0, 



'^A'^b^ 



2gr'^j 
0, 



The vacuum initial conditions are given by eq. l|10p : 

k 

'2\kV 



-00 



(28) 
(29) 
(30) 

(31) 
(32) 
(33) 
(34) 
(35) 

(36) 



The other components of the Wigner function have zero 
initial values. 

The system of equations in (P51) - ([55)) can be solved 
numerically for different color configurations. Thus the 
evolution of the color could be followed in detail, simi- 
larly to the color evolution in the Wong Yang-Mills equa- 
tions [3(1,^]. However, in the next section for the sake of 
simplicity we will consider assumption that the external 
field has a fixed direction in color space. 
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V. KINETIC EQUATION IN SU(2) WITH 
EXTERNAL FIELD OF FIXED DIRECTION IN 
COLOR SPACE 



The numerical solution of eqs. p4l) - p9)l requires initial 
conditions. In our case the appropriate vacuum initial 
conditions are the following: 



The equations ([28 l) -([35 l) provides a full description of 
the fermion production in a time-dependent SU(2) color 
field. Here we will consider the case of the external field, 
that has a fixed direction in color space: A" = A*n", 
where rt'^n" = I and 9tn" — 0, in general case, and = 
A^n"-, a = 1,2,3 in our calculation. Then the equations 
(P5)) - (|35p have the following particular solution: 

b'^ = h^n", d" = d^n". (37) 

Taking into account the zero initial conditions for 
and bQ,dQ, the equations ([32|) -(|35 |) become trivial and do 
not contribute to the set of equations (|28ll3ip . the last 
simplifies to give 



61(k,f = 



-oo 



— oo ) 
— oo) 



2k ' 

2k ' 
bl{k,t = 

d%k,t = 



-oo) = 
-oo) = 0. 



(51) 

(52) 

(53) 
(54) 



Now we have all parts to proceed numerically for any 
time-dependent field strength and determine the momen- 
tum distribution of the produced fermions. 



VI. NUMERICAL RESULTS 



d 
ok-i 



d_ 

d 



2k X d", 
2k X d^, 
2k X b^ 
2k X b^. 



(38) 
(39) 
(40) 
(41) 



Introducing the unit vector coUinear to the field direction, 
n = E*/|E*| = (0,0,1), and k^ = (fci,fc2,0), we can 
perform the following vector decomposition: 



b*l^ = bt-a + h 



d*l^ = d^l^n X ^ . 

kx 

Finally we obtain 6 equations, only: 



(42) 
(43) 



dtbl 
dtd'- 
dtb\ 



-E 



dk,^ 
3.9 d 



-E'^^d'' 



d_ 

'i dk^ 



dtd^ 



d 
oks 



These equations determine the time evolution of the 
wanted distribution function: 



//(k,t) 



fca 6| + ki_ 6] 



(50) 



In heavy ion collisions, one can assume three different 
types of time dependence for the color field to be formed: 
a) pulse-like field develops with a fast increase, which is 
followed by a fast fall in the field strength; b) formation 
of a constant field (Eq) is maintained after the fast in- 
crease in the initial time period; c) scaled decrease of the 
field strength appears, which is caused by particle pro- 
duction and/or transverse expansion, and the decrease 
is elongated in time much further than the pulse-like as- 
sumption. 

Figure [T] displays three sets for the time dependence 
of the external field I33l: 



e: 



pulse 



e^ 



it) =Eo 
tit) - 



[l-tanh^(V(J)] 
it) at t <0 



pulse 

Eo 



at t> 



(55) 
(56) 



= -2k3d' , 


(44) 




= 2kxd\ 


(45) 




= 2fc361 - 2fcj.6^, 


(46) 


-1— 

LU 


- -2fc3d^ 


(47) 




= 2k xd", 


(48) 




= 2/036! - 2fc_L63. 


(49) 






FIG. 1: The time dependence of external field E{t) in three 
physical scenarios: a) pulse {dotted line); b) constant field, Eq 
{dastied line); c) scaled decrease {solid line). 
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FIG. 2: The momentum dependence of the transverse, color 
singlet component b''j_{k±, k^) at t = for the Bjorken 

expanding scenario. 



FIG. 4: The momentum dependence of the transverse, color 
component b°j_{k±, fes) at t = for the Bjorken expand- 

ing scenario. 



(i+t/to)" at t > 

pl/2 



(57) 



In eq. ([55)) we choose S — 0.1/ i?g , which corresponds 
to RHIC energies. In eq. (|57|) the value k = 2/3 in- 
dicates a longitudinally scaled Bjorken expansion with 

1/2 

to = 0.01/ i?Q . In fact, the whole time dependence is 
scaled by Eq^^^. 

Before comparing the momentum distributions ob- 
tained numerically in the three different physical scenar- 
ios of eqs. ([55l) -(|57 |) we investigate our results on the quan- 
tities responsible for the momentum distribution. De- 



tailed numerical results are shown in the following five 
figures for the specific case of the Bjorken expansion. 

The distribution function ff{k±,k3) depends directly 
on b']^{k±,k3) and bl{k±,k^) as eq. ([SHI) shows. Figure 

and [3] display the magnitudes of these quantities in 2- 
dimensional plots. The fcs-symmetry of ^^^{k^, k^) can 
be seen clearly, as well as the asymmetric behaviour of 
6|(fcj^, fcs). The value of these functions are in 0(1). 

The distribution function ff{k±,k3) is color neutral, 
thus the color quantities 6^(fc_L, k^) and b'^{k±, k^) do not 
contribute directly. However, they play important role in 
the kinetic equation, see eqs. (|44p - P9|) . Their oscillating 
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FIG. 6: Longitudinal momentum (fea) spectra for fermions at 
kj_/Ef/^ — 0.5 and t = in the three physical scenarios 

(see Fig. 1 and text for explanation). 



FIG. 8: Transverse momentum spectra for fermions at 
k-i/El/^ = 0, 1, 2 in the time step t = for the Bjorken 

expansion scenario. 



behaviour can be seen on Figure 3] and O 

Now we return to the three physical scenarios with dif- 
ferent time dependence described in eqs. ([55)) - ([57|) and in- 

1 /2 

vestigate the fermion spectra at a large time, t = 2/Eq . 
Figure [6] indicates the behaviour of the longitudinal mo- 
mentum spectra at small transverse momentum value, 

1 /2 

where we choose k±/EQ = 0.5. Pulse- type time de- 
pendence leads to a narrow fcs-distribution {dotted line), 
which mimics a Landau-type hydrodynamical initial con- 
dition. The longitudinal spectra from constant field sce- 
nario {dashed line) leads to a flat distribution function in 




10-5 1 1 . 1 . 1 . 1 . 1 1 

-4 -2 2 4 

FIG. 7: Longitudinal momentum (fca) spectra for fermions at 
k±/Ey^ = 2.5 and t = 2/Eq^^ in the three physical scenarios 
(see Fig. 1 and text for explanation). 



k^. This result agrees well with a 1-dimensional, longitu- 
dinally invariant hydrodynamical initial condition, as we 
expect. Considering the scaled field scenario {solid line), 
since its time dependence is very similar to the pulse-type 
case (see Fig. [1} , then the similarity in the longitudinal 
spectra is well understandable. 

If we increase the transverse momentum and choose 
k^l Eq — 2.5, then Figure[7]displays the obtained longi- 
tudinal spectra: the ^3 dependence becomes very similar 
in the small ^3 region for the three different time evolu- 
tion. For pulse- type time dependence {dotted line) and 
for constant field scenario {dashed line) even the mag- 
nitude is very close to each other. Further numerical 
investigations are needed to understand this similarities. 

The interplay between transverse and longitudinal mo- 
mentum spectra is displayed on Figure [51 In the Bjorken 

expansion scenario of eq. (|57p we choose different longitu- 

1/2 

dinal momentum windows, namely k^/ E^ = 0, 1, 2 and 

1/2 

extract the transverse momentum spectra dX t = 2/i?Q . 
On Figure [5] one can see large deviation at small trans- 
verse momenta and similar spectra at high transverse 
momenta: hard fermion production is similar at different 
rapidities, but soft production (and any 'effective tem- 
perature') is very strongly rapidity dependent. 

Figure [9] displays the transverse momentum spectra for 
the three different physical scenario at momentum fc3 = 

1 /2 

and time t = 2/i?g' . Pulse- type time dependence leads 
to exponential spectra {dotted line), ff oc exp{kT /T) 

1 /2 

with slope value T — 1.54 • Eg . In the other two cases, 
we obtain non-exponential spectra generated by the long- 
lived field. Here the spectra from constant {dashed lines) 
and scaled {solid lines) fields are close, because the pro- 
duction and annihilation rates balance each other. Slight 
differences appear because of the fast fall of the scaled 
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10 20 30 40 50 



FIG. 9: The transverse momentum spectra for fermions at 
fcs = and t = in the three physical scenarios (see 

Fig. [T] and text for explanation). 



field immediately after t = 0. 

In our previous paper [33| we considered fermion pro- 
duction in U(l) classical field. Comparison between cur- 
rent SU(2) calculation and U(l) results can be made, 
if we fix both field strength at the same value. Tak- 
ing the value E^"^ / = 0.54 used in Ref. [Ill, we ob- 
tain transverse momentum spectra depicted in the Fig- 
ure [101 where particle production in the scaled field sce- 
nario was considered. We have found that the obtained 
transverse distribution functions for fermions are close to 
each other and have at almost the same shape. Figure [TT] 



displays the ratio of the two numerical results on a lin- 
ear scale. It clearly shows the appearance of a numerical 
value of /su(2)(fc-L)//u(i)(fcj-) = 0.75 for transverse mo- 
menta > As- This value may indicate the presence 
of scaling solutions of the investigated kinetic equations. 
The direct study of the kinetic equations in Section IV 
and V does not reveal scaling, but special cases can in- 
dicate such a behaviour. 



VII. SCALING SOLUTIONS 

In this section we consider specific initial conditions 
and corresponding solution for the kinetic equations 
eqs. (|44l) - (|49|) . At first we assume the singlet and the 
color multiplet components have the same initial values. 
This "symmetry" is generally violated by vacuum initial 
conditions demanding and (P to be zero in vacuum in 
contract to initial conditions for 6* and c?* (cf. ([5T|l - ((54l) ). 
Accepting this assumption, one may note that eqs. ([33])- 
(|49)) have the following special solution: b" = rib'^ and 
= rjd'' , where 77 is a numerical value to be defined in 
Appendix A. In this case the SU(2) kinetic equations for 
massless fermions can be shortened and rewritten into 
the following form (see Appendix A for details): 



V = VF(1 - 2/) - 2|k|M 
u = 2|k|i; . 



(58) 

(59) 
(60) 



One can recognize that formally these equations have 
been used earlier in the Abelian case to calculate massless 
fermion production (see eqs.(22)-(24) with zero fermion 
mass in Ref. 13311). 
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FIG. 10: Transverse momentum spectra of fermions from our 
calculation with scaled time evolution, Escaied{t), for SU(2) 
case (solid line), and from Ref. fs^ for U(l) (dash-dotted 
line). 



FIG. 11: The ratio of the numerical results from Figure lTOl on 
transverse momentum spectra for fermions in a linear scale: 
the SU(2) case (this calculation) is divided by the U(l) case 
from Ref. [H. 
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The region of our interest is the small longitudinal 
(fca —^ 0) and high transverse {k± —^ oo) momenta. It is 
reasonably to assume that for this region the distribution 
function is much smaller then unity, f{k± — > oo) <^ 1, 
demonstrating small particle yield in the high transverse 
momentum region. Thus the term of (1 — 2/) in cq. (j59|) 
can be simplified to unity. 

The eqs. (|58II60I) allow us to introduce new scaled vari- 
ables, namely W — Wrj, f — frj^, v = vrj and u — urj. 
Thus the above equations can be rewritten into the fol- 
lowing form: 



/ 

V 

u 



M^-2|k|u, 
2|k|i;. 



(61) 

(62) 
(63) 



In the U(l) case we have solved this set of equations with 
the scale parameters 77 = 1, and obtained the f{rj — 1) 
momentum distribution functions [ssl ]. Thus we know 
any scaled solution: /(ry) = /(ry = 1) • 77"^. In the case 
of SU(2) we have scale variable 77 = ±2/\/3 (see Ap- 
pendix A). Thus we can easily extract the wanted ratio: 
/su(2)//u(i) — 3/4. This number agrees with the numer- 
ically calculated value in the high transverse momentum 
region, as Figure [TT] displays. 

For small transverse momentum our analysis cannot 
be performed, because Pauli suppression factor (1 — 2/) 
differs from unity (see e.g. Figure O. The numerical 
calculations have shown a ratio of /su(2)//u(i) > 3/4 in 
this momentum region (see Figure [TT|) . 

The physical meaning of this ry-scaling is evident: 
number of particles with high transverse momentum is 
proportional to the second power of the field strength, 
namely {gE}'^. 

An additional scaling property has been used during 
the presentation of our numerical values. Namely in the 
eqs. (P5)) - (PT|) one can introduce a basic scaling in the 
variables: 



k 

Eo 
t 



A • fc , 

• So , 
A-i.<. 



(64) 
(65) 
(66) 



In this case the obtained distribution function will not 
change. 

Furthermore, in the case of external pulse field of 
eq. ((55)) . at t — > 00 the exponential transverse spectra 
of the stationary solution for the distribution function / 
can be characterized by an effective temperature as we 
have shown on Figure [HI This temperature will be scaled 
as the momentum, T — > A • T . 

In our previous paper [33| for the Abelian external 
field with scaled time evolution we have obtained numer- 
ical transverse spectra, which was very close to the per- 
turbative QCD results for high-fcr, namely scaling with 
log(fcT/As) • {kr/As)^ m. As Figures [TO] and [m dis- 
played, we obtain the same behaviour for the SU(2) case, 
although we can not prove this scaling analytically. 



VIII. CONCLUSION 

In this paper we investigated fermion production from 
a strong classical SU(2) field in a kinetic model. We 
derived the appropriate system of differential equations, 
which became simplified after introducing zero fermion 
mass. Following the physical ideas suggested in our ear- 
lier paper on the U(l) case (see. Ref. [3^]), we assumed 
three different time dependence for the longitudinal color 
field of fixed direction in color space and obtained three 
different types of longitudinal and transverse momentum 
spectra for the SU(2) fermions. We have found that the 
time dependence of the external field determines if the 
transverse momentum spectra are exponential or poly- 
nomial, similarly to the U(l) case. We solved the kinetic 
equation numerically and displayed the obtained results 
on longitudinal and transverse momentum distributions 
for the massless fermions. Furthermore, for the scaled 
time dependent external field in the transverse momen- 
tum spectra we obtained a constant ratio of 0.75 at high 
Pt between the recent SU(2) and earlier U(l) results. We 
could reproduce this scaling behaviour analytically in the 
high transverse momentum limit, together with the fac- 
tor of 0.75. Furthermore, our numerical results display 
the presence of a log(A:T/As) • {kr/Ag)^ scaling, similar 
to results obtained from perturbative QCD calculations. 



IX. APPENDIX A 

In this Appendix we display a specific solution of the 
SU(2) kinetic equation of 0, which leads to the form of 
equations known and solved already in the U(l) case (see 
Ref. nil). These equations allow us to recognize scaling 
properties in U(l) and SU(2), discussed in Section VII. 

Considering the form of SU(2) eqs. they can 

be solved assuming a very specific condition for the sin- 
glet and the color multiplet components [2^, namely 

fe^ = rib' , (67) 
d" = 7]d' . (68) 

Substituting this constraint into eqs. (|44p - ((49l) . a def- 
inite value can be obtained for parameter 77, namely 
77 = ±2/V3. 

In parallel, these equations become even simpler: 
d 



{dt + -E' 



(dt 
{dt 



9_ 

V 
9 
V 



E' 



dk3 

dks 
d 



)61 = ~2k3d% (69) 
)6^ 2fc_Ld', (70) 
)d' = 2k3bl - 2k±bl; (71) 



For further consideration it is appropriate to rewrite 
this set of equations in terms of the distribution function 
ff{k,t). Applying the differential operator 

77 9^3 



Dt^dt + ^E 



(72) 
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to the distribution function (|50|) and using (l69l) - ([7T|) we 
obtain 

Dtf = ^Wv, (73) 

Dtv = W{l~2f) ~2\k\u, (74) 
Dtu = 2\k\v; (75) 

where we introduced the functions W,v,u as 

W = (76) 

w ^ -2^". (78) 
Using method of characteristics [4l[ we finally get: 



/ - Iwv, 

V = VK(1 - 2/) - 2|fc|w, 
ii — 2\k\v; 

where ks becomes time-dependent 



(79) 

(80) 
(81) 

(82) 



One can recognize that the same equations have been 
used to calculate fermion production in the Abelian case 
with rj = 1 (see eqs.(22)-(24) with zero fermion mass 
in Ref. (33j). This surprising discovery indicates that in 
very specific cases Abelian-like equations can be derived 
from non-Abelian kinetic equations. 

However, there is a conflict between the constraint 
in eas. ((67)) - ([68|) and the physical vacuum solution in 
eq. ()10|) . which indicates a more complicated connection 
between the non-Abelian and Abelian cases. 
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